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Polyhedral vesicles with a large bending modulus of the membrane such as the gel phase lipid 
membrane were studied using a Brownian dynamics simulation. The vesicles exhibit various poly- 
hedral morphologies such as tetrahedron and cube shapes. We clarified two types of line defects on 
the edges of the polyhedrons: cracks of both monolayers at the spontaneous curvature of monolayer 
Co < 0, and a crack of the inner monolayer at Co > 0. Around the latter defect, the inner monolayer 
curves positively. Our results suggested that the polyhedral morphology is controlled by Co. 
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Amphiphilic molecules such as lipids and detergents 
form various structures such as micelles, cylindrical struc- 
tures, and bilayer membranes in aqueous solution [|l|. In 
particular, closed bilayer membranes, i.e., vesicles, are 
biologically important as model systems for the plasma 
membrane and intracellular compartments in living cells. 
Various morphological changes in the vesicles are under- 
stood via the Helfrich elastic model 0, § |. However, 
this model can not be applied to non-bilayer structures. 
For example, in an inverted hexagonal Hu phase, the 
hydrophobic interstice (void) space opens among three 
cylindrical monolayers. Recently, it is considered that 
this interstice space is filled by the tilt deformation of 
amphiphilic molecules as shown in Fig. |l|(a) The 
molecules tilt with respect to the monolayer surfaces 
around the junction of the three bilayer s. The monolayer 
surfaces are sharply bent at the junction. The effects of 
tilt deformation are also studied regarding the fusion in- 
termediates of fluid phase membranes 

On the other hand, polyhedral-shaped vesicles of fim 
scale size were observed in the gel phase; a triangular 
pyramid or prism-shaped vesicle of a mono-component 
lipid [D , and an icosahedral vesicle of mixtures of cationic 
and anionic surfactants 0. The membranes are flat on 
the faces of polyhedrons and are sharply bent at the 
edges. Because the bending modulus are very large in 
the gel phase, the polyhedral vesicles would be stabler 
than the spherical vesicles. The free-energy loss of the 
defects at the edges would be less than the loss of the 
equal bending of membranes on the sphere. However, 
the defect structure at the edges is unresolved. The in- 
formation of the edge structure is significant to control 
the morphology of the polyhedral vesicles. These vesicles 
are expected to be of practical value for drug delivery. 

To clarify the edge structure, theories or simulations 
with molecular resolution are needed. Since molecu- 
lar dynamics simulations with atomic resolution have 
been applied only for the ~ 10ns dynamics of 1000 
lipid molecules due to the restrictions of computational 
time |8|, coarse-grained molecular simulations ||, |l^, 
pj] , |l2, |l3| have been applied. We studied the fusion 
and fission dynamics of vesicles using Brownian dynam- 




FIG. 1: Three types of Line defects, (a) Tilt deformation in 
the inverted hexagonal phase, (b) Cracks of both monolayers. 
Hydrophobic segments are partially exposed, (c) Crack of in- 
ner monolayer. Amphiphilic molecules in the inner monolayer 
tilt with respect to the boundary surfaces of two monolayers. 



ics [|ll[. The self-assembly into vesicles is simulated by 
our model a lattice Monte Carlo method Q, and 
dissipative particle dynamics Q. However, these simu- 
lated vesicles were flexible, and no polyhedral vesicle has 
been obtained. 

In the present paper, we developed our previous model 
to control the bending modulus of monolayers by the 
addition of the curvature potential of a monolayer. We 
simulated the polyhedral vesicles of rigid membranes, and 
obtained two types of defects at the edges as shown in 
Figs, [^(b) and (c). The morphology of the polyhedral 
vesicles and the defect type depend on the spontaneous 
curvature of the monolayer Co- 

An amphiphilic molecule is modeled as rigid rods con- 
sisting of one hydrophilic segment (j = 1) and two hy- 
drophobic segments {j = 2,3), which are separated by 
a fixed distance a. Solvent molecules are not explic- 
itly taken into account, and "hydrophobic" interaction 
is mimicked by the multibody local density potential of 
the hydrophobic segments. As details of the basic model 
were described in our previous papers [|l^, |ll|, here we 
briefly explain the model. The motion of the jth segment 
of the ith molecule follows the underdamped Langevin 
equation. Amphiphilic molecules {i = 1,..,A^) interact 
via a repulsive soft-core potential t/rep, an attractive "hy- 
drophobic" potential /7hp, and a curvature potential Ucv- 
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FIG. 2: Sliced snapshots of vesicles at temperature k^T je — 0.2 and the number of molecules N — 1000. (a) Disk shaped 
(dihedral) vesicle at the spontaneous curvature of monolayer Cocr = —0.11. (b) Triangular- pyramid shaped (tetrahedral) vesicle 
at Cocr = 0.058. (c) Pentagonal-prism shaped (septihedral) vesicle at Cqct = 0.23. Gray spheres and white cylinders represent 
hydrophilic and hydrophobic segments of amphiphilic molecules, respectively. 



where pi^j is the local density of the hydrophobic 
segments for the jth segment of the zth molecule. 
Both segments have the same soft radius, ?7rep(^)/^ = 
exp{— 20(r/<j — 1)}. The "hydrophobic" interaction is 
mimicked by the potential /7hp(p). 

= XI ^(I^M -ri'j'l), (2) 

where hir) = r— — — — — — ^ -. 

^ ^ exp{20(r/cr- 1.9)} + 1 

Pi J is the number of hydrophobic segments in the sphere 
whose radius is approximately 1.9<j. 

r -o.5p {p<P*- 1) 

Uhp{p)/e = { 0.25(p - p*)^ -c {p*-l<p<p*) , 
{ -c {p* < P) 

(3) 

where c is given byc = 0.5p* — 0.25. We used the values 
p* = 10 and 14 at j = 2 and 3, respectively. 

To give the bending modulus n and the spontaneous 
curvature Co of the monolayer membranes, we use the 
potential /7c v of the orientational difference of neighbor- 
ing molecules. 

Ucv = X0.54Xr,,,0(ui - u,. - C^V)', (4) 

where the vector is the unit orientational vector of 
the zth molecule, and Yi^i' {vi^i') is the unit vector (dis- 
tance) between the ith and i'Xh molecules: = (r^^i — 
ri,3)/|ri,i - r^^sl and Yi^i> = (r^ - ri/)/|ri - r^/ 1, where 
is the center of mass of the ith molecule. At Cq = 0, 
this potential is similar to the bending elastic potential 
used in the tethered membrane models [Q. When the 



orientational vectors are equal to the normal vectors 
of the monolayers with no tilt deformation, 

Ucv/e = j 0.5a>:cv[(Ci+C2-Co)'-2CiC2 + Co']c^A (5) 

in the continuum limit, where Ci and C2 are the two 
principal curvatures of a monolayer. The spontaneous 
curvature Co equals to Cgcr/rnb, where fnb is the mean 
distance between neighboring molecules and fnb = 1.5cr. 
On the assumption of the hexagonal packing of molecules 
in the monolayers, we obtain = We used 

= 3^ to represent the rigid membrane, and /^cv — 
In previous papers, we estimated the bending modulus 
hzo/e :^ 0.5 ( half of the bending modulus of bilayers) 
at hi'^^ = [|ll[. Since hicv is ten times larger than /^q, 
the bending elasticity is mainly given by Ucv^ and the 
bending modulus of the monolayer 

We mainly used the number of molecules N = 1000 
and k^T / e = 0.2, where /cb is the Boltzmann constant 
and T is the temperature. Amphiphilic molecules spon- 
taneously form vesicles in a fluid phase at k^T/e = 0.2 
and K,'^^ = 0. The unit length a corresponds to ~ Inm. 
The unit time step tq = C^cP' js corresponds to ~ Ins 
estimated from the lateral diffusion constant of phospho- 
lipid at 30°C, - 10-^cmVs H, where C is the friction 
constant of the segments of molecules. 

Vesicles exhibit various polyhedral morphologies at 
/^^^ = Ze. The number of faces nf of polyhedron increases 
as Co increases. Figure || shows examples of the polyhe- 
dral vesicles. The edges of the polyhedrons are formed 
by the line defects [Figs, [^(b) and (c)]. The molecules 
at the line defects are shown using the number of neigh- 
boring molecules 

nfb in Fig. I The inner monolayers 
are divided into nf faces by the defects. At Cq > 0, the 
cracks of the inner monolayer [Fig. |l|(c)] occur on the 
edges, and the outer monolayer consists of one curved 
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FIG. 3: Line defects (edges) of the tetrahedral vesicle in 
Fig. ^b). The snapshot is viewed from the same viewpoint. 
The molecules with number of neighboring molecules nf^ < 
4.5 in the inner monolayer are shown. The number nf^ is 
defined as n^^ = 
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FIG. 4: Spontaneous curvature Co dependence of the mean 
number of faces < rif > of polyhedrons at = 1000 and 
/cbT/s = 0.2. Annealing: Vesicles are annealed from k^T je — 
0.5 to 0.2. Quenching: Vesicles are simulated starting with 
spherical vesicles at — 0. a^cv increasing: the coefficient 
/^cv slowly increases from n'c^/e = to 3 starting with the 
spherical vesicles at a^cv = 0. 

face [Figs. ||(b) and (c)]. At Co < 0, the outer monolayer 
also exhibits cracks [Fig. ||(b)] on the edges. When the 
deformation of the left or right side of Fig. ^(a) is formed 
on the entire circular-line defects, the outer monolayer of 
the disk-shaped vesicles is divided into three faces (two 
disks and one cylinder) or two faces, respectively. 

Figure ^ shows the mean number of faces < nf > at 
k^T/e = 0.2 using three methods. The results of anneal- 
ing are the closest to the equilibrium values. Through 
the other methods, however, vesicles are often trapped 
in met ast able states. At Co > 0, vesicles with larger 
(smaller) nf values are obtained through quenching {hCcv 
increasing) than through annealing. Thus, one should 
use the annealing method to obtain regular polyhedrons. 

Flip-flop motion, which is the transverse motion be- 
tween the inner and outer monolayers, frequently oc- 
curs at k^T / £ = 0.5. The number of molecules in 
the inner monolayer decreases with an increase in nf at 




FIG. 5: Spontaneous curvature Co dependence of (a) the 
mean curvature for amphiphilic molecules, < Cam > (b) the 
mean curvature for monolayer surface < Csf > and the differ- 
ence < Cdif >=< Cam -Csf > at = 1000, kBT/s = 0.2, and 
7in = 0.328(=b0.003). The superscripts 'in' and 'out' represent 
inner and outer monolayers, respectively. Circles: rugby-ball 
shaped trihedron. Squares: tetrahedron. Triangles: triangu- 
lar prism. Diamonds: cube. 



k^T/e = 0.5, since tilting molecules on the line defects 
share a larger area. The ratios 7in of molecules in the 
inner monolayer are 0.31(±0.01) and 0.292(±0.003) at 
Cocr = —0.11 (nf = 2) and Cocr = 0.23 (nf = 7.3), respec- 
tively. On the other hand, flip-flop motion rarely occurs, 
and the ratio 7in is fixed at 7in = 0.328(±0.003) at the 
fixed temperature k^T / e — 0.2. In typical experimental 
conditions, flip-flop motion is very slow, and the half-life 
is more than several hours even in the fluid phase [Q. 
Thus experimentally, the ratio 7in of the polyhedral vesi- 
cles should not reach an equilibrium value as well as the 
simulation with the fixed temperature /ceT/e = 0.2. 

In the obtained polyhedrons, three edges are connected 
at any vertex. The connections of more edges are un- 
stable and are not formed. The number of edges Uq 
then equals to 1.5 times the number of vertices nv, be- 
cause each edge contacts two vertices. We derived Uq = 
3(nf — 2) and — 2{nf — 2) from this relation and Eu- 
ler's formula for a convex polyhedron (nf + nv — ne = 2). 
At nf > 6, multiple types of polyhedrons with the same 
number of faces exist. However, we obtained only one or 
two types of polyhedrons; i.e., cube (4^) at nf = 6; pen- 
tagonal prism (4^5^) and (3\4^,5^) at nf = 7; {4^,5^) 
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at rif = 8; (4^5^) and (4^,5'^,6^) at nf = 9, where (p^) 
represents a polyhedron with q p-gons. Thus more sym- 
metric polyhedrons are more frequently formed. 

Figure ^(a) shows the Co dependence of the curva- 
tures of four polyhedrons at Co > 0, where Cam is the 
curvature (splay) for amphiphilic molecules, < Cam > 
^nb = - Ui^)ri^i^h{ri^i'))/{^h{ri^i^)). The mean 

curvature < C°^* > of the outer monolayer is almost in- 
dependent of Co and nf. On the other hand, < C^^ > of 
the inner monolayer increases with an increase in Co and 
nf, since the inner monolayers curve positively around 
the crack of the inner monolayers. The mean curvature 

< Cam > foi" a larger nf is closer to Co, though more 
hydrophilic segments contact the hydrophobic segments 
at the line defects. The curvature Cam does not coincide 
with the curvature of the monolayers, since molecules can 
tilt with respect to the monolayer surface. To clarify this 
tilt, we estimated the curvature Cgf for the monolayer sur- 
face, < Csf > fnb = {Y,{ni-nif)Ti^ifh{ri^i^))/{Y,h{ri^i')), 
where is the normal vector of the monolayer sur- 
face at r^. We defined as the vector minimizing 

= X] when nf^ > 2.5. This mini- 

mizing vector is the eigenvector with the smallest eigen- 
value of the moment tensor of inertia of the neighboring 
molecules. Figure ||(b) shows the curvatures < C^f > of 
the inner monolayer surface and difference < C^^f > be- 
tween the two curvatures. The inner monolayer surface 
(molecules in the inner monolayer) tilt with respect to 
the boundary surfaces of the two monolayers (the inner 
monolayer). Both tilts increase Cam- Since the molecular 
tilt in the inner monolayer < C^^f > is almost indepen- 
dent of nf, the length of the line defects only changes the 
curvature < Cgf > of the inner monolayer surface. Thus 
the polyhedral morphology at equilibrium should be de- 
termined by the effects of line defects on < C^f > and 
the hydrophobic interaction. 

We found that the spontaneous curvature Co can be 
estimated using this tilt deformation with respect to the 
monolayer surface. At Cgf < Co (Cgf > Co), molecules 
tilt to reduce |Cam — Co| and Cdif shows a positive (neg- 
ative) value. Thus, we estimate Cocr ^ 0.04 at hz'^^ = 
and kBT/s = 0.2: < Cdif >cr = 0.0041 (±0.0002) for flat 
membranes; < C^J}* > a = 0.0(±0.00007) and < C^^* > 
a = 0.0429(±0.0005) for tube-shaped vesicles with a di- 
ameter of 30cr; and < C^f^ > a = -0.0056(±0.0003) and 

< Cg7* > cr = 0.0975(±0.0003) for spherical vesicles with 
a diameter of 20cr. 

The line defects may be interpreted using the correc- 
tion terms of the Helfrich model, the local minimum at a 
large Ci + C2 . The morphology of the polyhedral vesicles 
may then be obtained from the Euler-Lagrange differen- 
tial equation. Similar deformations to the crack of the 
inner monolayer are seen in our daily experience. When 
a rubber hose is strongly bent, a side of the hose becomes 
hollow, and the other side smoothly bends. 

The morphology of the polyhedral vesicles depends on 



the number N and properties of molecules. At N = 2000 
and Cocr = 0.23, vesicles exhibit more complex morpholo- 
gies with concave edges, where the membrane bends out- 
side with a crack of the outer monolayer. The hexagonal 
packing of molecules should stabilize the triangular and 
hexagonal faces. In some multi-component vesicles, the 
phase separation occurs at the edges or vertices of the 
polyhedrons. Dubois et al reported that the segregated 
anionic surfactants form pores at the vertices Q. Var- 
ious polyhedral vesicles are likely to be experimentally 
observed under the control of Co and other conditions. 
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